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1. Introduction 

In a renormalizable field theory, given enough time and computer power, 
one may calculate an observable up to any order n in the perturbative expan- 
sion, study the convergence of the series and make a precision comparison 
with experimental results. In the absence of complete higher order calcu- 
lations, a first estimate of the convergence of the perturbative series may 
come from the so-called leading logarithms (LL). These are terms of the 
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form a n log ra /i 2 , with [i the renormalization scale, which appear in the n-th 
order corrections of any observables expression upon renormalization. The 
coefficients of these LL may be calculated using renormalization group meth- 
ods. 

In a non-renormalizable field theory the situation is more convoluted, 
since new terms appear in the Lagrangian at each order in the expansion. 
However, as Li and Pagels pointed out [lj to one-loop, the n-loop-order con- 
tributions will still contain a c n log n (// 2 /M 2 ) term. Consider for example 
the expansion of the pion mass in the quark masses in Chiral Perturbation 
Theory 0]: 



Ml = M 



M2 Mi M 2 or, \ 



(1) 



where l\ is a renormalized second order Lagrangian coupling and M and F 
are the lowest order Lagrangian parameters. Depending upon the /i 2 /M 2 
ratio, where the parameter M 2 is the lowest-order mass and \x is the renor- 
malization scale, the term with the logarithm can be the largest part of the 
correction. In general the LL, now depending on a typical scale M. of the pro- 
cess one looks at, may turn out to be a substantial fraction of the n-th-order 
correction. For many observables indeed the LL are the main contribution 
due to the enhancement of log(/z 2 /M 2 ) compared to the other contributions 
at the same order. This is true in particular for the 5-wave scattering 
length a[] 0, 0| in Chiral Perturbation Theory (ChPT). 

Many corrections of this type at one-loop have been calculated long ago, 
see the review Q and references therein, mostly in the framework of Current 
Algebra, and in ChPT @, 0, @|. 

Weinberg pointed out in 1979 Q that, because of renormalization group 
equations (RGE), the two-loop LL coefficient appearing in tttt scattering 
amplitude could be calculated using simple one-loop diagrams. This method 
was later used in Q for two loop LL corrections to scattering lengths and 
slopes in mr scattering and in [8j for the general three-flavour meson sector. 
Nowadays the extension to the full two-loop expressions for mesons is known 
for most observables [9(. 

Weinberg's renormalization argument was extended by Biichler and Colan- 



gelo [10[ to all orders and to a generic non-renormalizable theory. They 
showed that the leading logarithms at any loop-order can be calculated us- 
ing one-loop diagrams. They also showed that the coefficient of the leading 



2 



logarithm only depends upon the constants appearing in the lowest-order La- 
grangian. In principle this coefficient at n-loop-order could have depended 
on all of the coupling constants in the Lagrangians C m with m < n. 

However, the problem remains that as n grows, the number of terms and 
counterterms in the Lagrangian grows very rapidly, and the renormalization 
group equations (RGE) become more involved. This renders the calculation 
of LL beyond the first few orders a Herculean task0. 

The alternative to performing these long calculations is to extract the LL 
series from a renormalizable theory in hopes that it will reproduce the LLs 



of the non-renormalizable theory. The authors of [12[ applied this approach 
to the renormalizable linear sigma model and were able to resum the entire 
LL series by exploiting recursion relations. They however found recursive 
relations were not possible in the non-renormalizable non-linear sigma case 
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In the massless case, a solution to managing the terms in the Lagrangian 
was found since the number of meson legs on the vertices one needs to con 



sider remains limited [13J. This was used for meson- meson scattering [13 



and the scalar and vector form factor 14| in the massless O(N) model. This 



method works to arbitrarily hig h order and agrees with the known large 



N [15| results to all orders [13j, ll4(. In the massless case, one may also 



use kinematic methods to extract the nonanalytic dependence on kinematic 



quantities. These have for instance been used to derive the form factors [16 



up to five loop-order and to arbitrarily high order for both the form factors 



and the meson-meson scattering amplitude [171 ] . These methods essentially 
solve the leading logarithm problem in the massless case for the most useful 
observables. 

In the massive case (which includes ChPT), however, tadpole diagrams no 
longer vanish and one needs to consider terms with an increasing numbers of 
meson legs. E.g. for the meson mass one needs to calculate one-loop diagrams 
with 2n meson legs in order to get the n-loop-order LL. In our earlier work 



18| we showed that one does not need to explicitly construct the higher order 
Lagrangians in a symmetric form, nor does one need a minimal Lagrangian 
at each order. The LL series only requires order by order a complete enough 



Lagrangian to describe the observable at hand [18|. This means that one 



may let the algorithm itself generate all the necessary terms in the higher 



^ee for example the two- loop leading logarithm in the non leptonic sector 11 1 
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order Lagrangians. 

We applied this method to obtain the leading logarithm to five-loop- 
order for the meson mass in the nonlinear 0(iV)-model. One set of results 
of this paper is to extend the calculation to the decay constant and vacuum 
expectation value also to five-loop-order. A similar amount of work allows 
to obtain the leading logarithms for meson-meson scattering and the scalar 
and vector form factors to four-loops and we present results for this as well. 
For N = 3 the massive non linear O(N) model is equivalent to two-flavour 
mesonic ChPT Lagrangian at lowest order in the sense that 0(4)/ 0(3) is 
isomorphic to 577(2) l X SU (2)r/SU (2)y ■ We have thus obtained the leading 
logarithms also for this physically interesting case to rather high loop-order. 



In [18j we also extended the large N limit to the massive case and applied 
it for the meson mass. In this paper we extend those methods to many 
more observables and we find for all the cases considered simple expressions 
in terms of t he p hysical quantities. However, as already observed in the 



massless case [13, |lj] and for the mass 18j , we again find that the large N 
result does not give a good approximation to the coefficient of the leading 
logarithm for general N and in particular not for N = 3. 

We briefly summarize the methods for calculating leading logarithms of 
1 181 ] in Sect. [2] and introduce the massive 0(N) nonlinear sigma model 



including external fields in Sect. |3j Here we also define all the physical quan- 
tities we calculate. The next section discusses the large TV-limit. We briefly 



recall the results of [18] and extend the method to the other observables. 
The leading logarithms for general N are discussed in Sect. |5j There we also 
discuss the convergence of the various observables and compare some ways 
of expanding. The main results and conclusions are summarized in the last 
section. 

One note about the cross-checks on our results. All calculations were 
performed using four different parametrizations for the fields. This means 
that for every parametrization the form of the Lagrangian and the couplings 
are different. These four different Lagrangians were fed into the same form 



FORM code [19|. The fact that the output for all observables came out the 



same regardless of the parametrization is a very good sanity check. 
2. Counter terms and Leading Logarithms 



In this section we present the results of [10] and [18j. We show how to 



calculate the LL coefficients and their connection to the counter terms used 
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to renormalize the observable. 

In order to calculate the matrix elements for the observables, consider the 
generating functional 

W\j\ = e iZ[j] /h = [ V^e iS ^ n . (2) 



where j is the classical source which allows one to extract all of the Green 
functions. The action S may be expanded around the classical action, S = 
En=o hn S {n) - In practice one expands the Lagrangian into a sum of growing 
h order (bare) Lagrangians 

oo 

£ = J2 h n C {n) . (3) 

n=0 

The crucial difference between renormalizable and non-renormalizable theo- 
ries is the number of terms appearing in each = Ef"4 n) Oj n) . In the 
former case the terms in all & n > are of the same form as in £(°). In the later 
case, new terms appear at each order n. When calculating matrix elements 
beyond tree level, the loop corrections lead to divergences, which must be re- 
absorbed into a redefinition of the coupling constants. While for a renormal- 
izable theory it suffices to reabsorb them into lowest order Lagrangian since 
they are all of the same form, in a non-renormalizable theory a divergence is 
absorbed by the higher order coupling constants. E.g. the divergence that 
appeared in the calculation of ([I]) was absorbed into £3. Alternatively, one 
may say that the renormalization consists in adding to the Lagrangian order 
by order operators with diverging coefficients to cancel divergences, i.e. 
counterterms c^O^ /e k . 



£( n ) _ ^•(n),ren _|_ £(n),di\ 

/'" 



1 



k=i e 



(4) 



where we have assumed one works in dimensional regularizatiorfl in which 
divergences appear as poles l/e fc , e = 4 — d. We have shown in [l8| that set 
operators need not be minimal or even complete for our purposes. 

Consider an observable and the one particle irreducible (1PI) diagrams 
that may contribute at order n. That one only needs to consider 1PI diagrams 



2 Thc pre-factor ensures that all Lagrangians have the same dimension d. 
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was proven in 10] . Let be the n-th order contribution from I loops. 



consists of a finite part and (possibly) several different poles 

l j n 

LI = + E ~f- (5) 

fc=i e 

In L™ k /e k we include only those divergences coming from the loop integration 
and not those coming from the diverging 0{ vertices in the loops. Each 
loop may contain different c^Oi vertices. We indicate with {c}" the product 

c hki c i2k2 ' ' ' c iTkr giving an n-th order, £-\oop contribution. 

The recursive equations follow from the requirement that the divergences 
must cancel. The contribution at order h may be written as 



+ -A» (too) + ^o- (6) 



where L\ x = L\\ ({c} ,). Expanding // 6 = 1 — elog^u + • • •, one finds that to 
cancel the 1/e we need 

Ll ({c}\) = —L\\ . (7) 



This determines all the needed in terms of the lowest order parameters. 
From ()6]) we also get the explicit dependence on \i 

-Lj ({c}J)log^=Lj 1 bg/i (8) 

where the equality follows from ([7]). To summarize, the counterterm c[i is 
adjusted so that it cancels the divergence coming from the loop, L\ x , this in 
turn determines the \x dependence and hence the coefficient of the LL. At 
second order the cancellation of the 1/e 2 and log(//)/e pieces allow to obtain 
the leading divergence from the one-loop part -^ii({c}J) and expanding 
/i e one finds that the coefficient of the LL log 2 (/i), is L\ 2 . See [l8[ for a more 
detailed discussion. 

These results may be generalized. At order n one may write 



1 



(9) 
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Requiring that the coefficients of 1/e™, log /// e n 1 , log 2 jj/e n 2 , . . . cancel leads 
to a set of n equations, the solution of which is given by [l8[ 

K-i,n-i {{C}f) = ( ) Lin- (10) 



In particular, 
and 



L n n[{c} n n Z\ =(-!)" nL nn (11) 



Ln ({c} n n Z\) = (-n)L n Q0 (12) 
The coefficient of the leading logarithm is given by 

^(log/if. (13) 



Eq. (TTTj) is solved recursively. First one calculates the one loop counterterm. 
With this one, using f JTTj) . one may calculate Z/ 2 2 , the coefficient of the second 
order LL. This again fixes the counterterm, which can be inserted back 
into Eq. (iTTj) . and so on. One only needs to insure that all the c^ n }\ for all 
the 0} n-1 ^ appearing in the calculation are determined. 

3. Massive nonlinear 0(N + 1)/0(N) sigma model 

The 0(N + 1)/0(N) nonlinear sigma model, including external sources^), 
is described by the Lagrangian 

C na = ^-D^ T D^ + F 2 x T $ , (14) 



3 If one wishes to study a given current one adds an extra classical source field to 
the generating functional which couples to that current. Thus the generating functional 
becomes 



The matrix elements involving J M can be obtained by functional derivation with respect 



to 
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where $ is a real N + 1 vector, $ T = $ x . . . $ N \ which transforms 
as the fundamental representation of 0(N + 1) and satisfies the constraint 
$ T $ = 1. The covariant derivative is given by 

D^ a = d^ a + v ab $ b - a^$° . (15) 

The vector sources satisfy t>|f = — v h ^ and correspond to the unbroken group 
generators while the axial a° sources correspond to the broken ones. Indices 
of the type a,b, . . . run over 1, . . . , iV in the remainder. The mass term x T & 
contains the scalar, s°, and pseudo-scalar, p a external fields as well as the 
explicit symmetry breaking term M 2 . 

T _ ( /od„0 , 7.^2^ 1 N 



X 1 = { (2Bs u + M 2 ) p L ... p N ) . (16) 

The term proportional to M 2 breaks the 0(N + 1) symmetry explicitly to 
the O(N), whereas the vacuum condensate 

($ T > = (1 ... 0) , (17) 

breaks it spontaneously to the same O(N). 

This particular model corresponds to lowest order two-flavour ChPT for 
N = 3 P, |20|. It is also used as a model for a str ongly i nteracting Higgs 



sectors in several beyond Standard models, see e.g. |21|, |22 . 

The terminology for the external sources or fields is taken from two-flavour 
ChPT. The vector currents for iV = 3 are given by v ab = e cah v c with e cba the 
Levi-Civita tensor. The electromagnetic current at lowest order is associated 
to v 3 . The external scalar source s° is instead associated to the QCD current 
— uu — dd as in (2)]. 

We write $ in terms of a real iV-componentQ vector (p, which transforms 
linearly under the unbroken part of the symmetry group, O(N). In the calcu- 
lations, for simplicity, we will refer to one particular parametrization, called 
$1 below. We have, however, made use of four different parametrizations in 



4 We refer to these as a flavour components. 
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order to check the validity of our results. These are 



$1 



1 - 



F 2 



$9 



1 



4>i $ 



1 V <P 

2 F 2 



1 T <p 
4 F 2 F 



F 

COS 



U) 

F 2 



\ 



sin 



F 2 



/ 



$i is the parametrization used in [2], $ 2 a simple variation, $ 3 is such 
that the explicit symmetry breaking term in (I14p only gives a mass term to 
the field but no vertices. $4 is the parametrization one ends up with if 
using the general prescription of 23]. 

The physical mass of the meson (squared), Mp hys , we already calculated 
in [isj. The meson decay constant, -F p h ys , is defined by the matrix element 
of the axial current j b a 



be 



(0|il^(P))=^physP^ 



The lowest order is Fphys = F. 

The vacuum expectation value (VEV) is defined by 



(18) 



V T 



phys 



(-J° s o) N = 3 (Uu + dd) 



(19) 



In the second equation we have written out its meaning in terms of quarks. 
Its lowest order value is Vphys = Vq = —2BF 2 . 

The scalar form factor is defined as the matrix element of the scalar 
current between two mesons 



IP/J 



i s °o| 



*»(p0) = F 8 [(p f - 



(20) 



The value at zero momentum transfer can be derived from the meson mass 
via the Feynman-Hellmann theorem as 



F s (0) = 2B- 



dM- 



phys 



DM 2 



(21) 



The vector form factor is defined similarly as 
(0°(P/)b'i-jt M l0 6 (Pi)} = (6 ac 5 db -5 ad 6 bc )z(p f + Pi yF 



v 



vP/ - PiJ 



(22) 
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where we have exploited the antisymmetry of v® b . The vector currents j v for 
the N = 3 are given by j ab = jyS cab with e cab the Levi-Civita tensor. The 
electromagnetic current in this case to the lowest order is given by j Vfl = 
{vrfpU — d r y fl d)/2. The singlet part does not appear to lowest order. The 
value of the vector form factor at (p/ — p^) 2 = is always 1 because of the 
conserved O(N) symmetry. 

In addition to the form factors we also define the radii and curvatures 
with t = (p/ — p^ 2 and expanding around t — 0: 

F s (t) = F s (0)(l + ^(r 2 } s t + c s t 2 + ---y 

F v (t) = l + l(r 2 ) v t + c v t 2 + ■■■ . (23) 
o 

The final process we discuss is meson-meson scattering. The general 
amplitude is 

(0 a (Pa)0 b (p b )|0 c (p c )0 d (Pd)) = 5 ab 5 cd A(s,t,u) + 5 ac 5 bd A(t, Ul s) 

+5 ad 5 bc A(u, s, t) , (24) 

with 

S = (Pa + Pb) 2 , t = (p - p c ) 2 U = (p a p d ) 2 , (25) 

satisfying s + t + u = 4M 2 hys . A(s, t, u) is symmetric in t and u. The proof 
of (j24p for = 3 is done using crossing and isospin symmetry 24j, but may 
be generalized to the O(N) symmetry case. There are three channels, the 
singlet, antisymmetric and symmetric combination, named / = 0, 1 and 2 
(after isospin conservation). The amplitude in these three channels is given 
by 

T°(s,t) = NA(s,t,u)+A(t,u,s)+A(u,s,t) 
T^Sjt) = A(t,u,s) - A(u,s,t) 

T 2 (s,t) = A(t,u,s)+A(u,s,t). (26) 
These amplitudes are projected onto the partial waves 

T l = 7TT- t d(cos6)P e (co S 6)T I (s,t) , (27) 

with 9 the scattering angle and Pi the Legendre polynomials. Near threshold 
the partial waves can be expanded in terms of scattering lengths a\ and slopes 
b\. 

n [T/1 = q 2e U\ + b\q 2 + ■••), (28) 
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where q 2 = | — 4M 2 hys ) . The scattering lengths and slopes are normally 
given in units of powers of M p h ys . 

For all the quantities defined here, the N = 3 results with a complete 
0(p 6 ) Lagrangian are known up to two loops, thus including the leading 
logarithms. This is an additional check on our calculation. 

The results can be expressed in terms of the lowest order parameters, 
expanding in powers of 

M 2 jj 2 

log tW , (29) 



16tt 2 F 2 to M 2 

or in terms of the physical mass and physical decay constant using 

Lphys = le^F 2 log ~W~ • (30) 

1D7T -Tphys iW phys 

In both cases we chose the mass scale in the logarithm to be the corresponding 



mass. In [18j we also used an expansion in 



H, s -^ lo g^' (31) 

We use this hybrid form in this paper only for one figure. 

4. Large iV limit 

When we consider the limit of large N we have to decide how the constants 
in the Lagrangian (fT4"j) depend on N. The first term can be made linear 
(extensive) in N by assuming F 2 cc N since $ T $ = 1. The second term is 
then also linear in N if we set M 2 and B to be independent of N. 



In the linear O(N) model it is well known [15|, |2lj that the leading con- 
tribution in N comes from diagrams that contain only non-overlapping loops 
and in which each factor 1/N coming from a new interaction is canceled by 
the factor iV introduced by summing over the N internal mesons in each 



bubble. In 18] we showed how this is also true in general for the non linear 



case which has vertices with any number of meson fields. The proof in [18 
remains valid when vertices with external fields are included. The vertices 
with more mesons legs are suppressed by powers of 1/F as was the case for 
the purely mesonic vertices. The suppression factor compared to the lowest 
order of 1/ F 2L with L the number of loops remains thus valid as well. 



11 



Figure 1: A typical diagram which contributes at leading order in N. Note that vertices 
can have many different loops attached since the Lagrangians contain vertices with any 
number of fields. The flavour-loops coincide with momentum loops, i.e. the visible bubbles. 




Figure 2: The graphical representation of the equation that generates all the cactus dia- 
grams for the propagator. A thick line indicates the full propagator, a thin line the lowest 
order one. 



In Fig. [T] we show such a typical diagram contributing to the meson 
self energy. All of these diagrams can be resummed by exploiting recursive 
methods. 

Consider for instance the inverse of the full propagator, it is given by the 
inverse of the lowest order propagator and the sum of all the 1PI diagrams 
with two external legs. By starting out with the lowest order propagator 
on the right hand side (rhs) and then reinserting the solution recursively we 



generate all cactus diagrams. In [18] we used this method to show that the 



full propagator in parametrization 1 in the large N limit is 

*A full (q 2 ) = 1 . (32) 

H JW phys 
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-Q/\/v • 




Figure 3: The graphical representation of the equation that generates all the one-particle- 
irreducible cactus diagrams for the decay constant. A thick line indicates the full propa- 
gator. The photon line indicates the insertion of the axial current a^. 



Note that, as shown in [18( and below, in this parametrization in this limit 
there is no wavefunction renormalization. Let us note once more that in the 
following we derive the results in the first parametrization but the others can 
also be used and give the same results. 



This method is similar to the gap equation used in e.g. NJL models [25 
but we have here an infinite number of terms on the right hand side. Similar 
resummations may be used for the other observables as shown below. 



For completeness we quote the result for the physical mass [18j], it is the 
solution of 

M 2 = M p 2 hys ^l + ^(M p 2 hys ) , (33) 

with 

= is? (34) 

Note that in non renormalizable field theories the large N limit also de- 
pends upon how depend on N. The result (|33l) and those derived below 
assume that the finite part of any higher order coefficient vanishes, c-q >0 ^ = 0, 



--iO 

see 26( for a discussion. The formulas still give the leading logarithms even 
for non zero higher order coefficients as long as the Lagrangian remains at 
most linear in N. 



4-1. Decay constant 

From the Lagrangian (I14p we can extract the vertices involving the axial 
current j£. These are given by F 2 (9 /i $°a°$ a — <9 M $ a a^$°). In the first 
representation $ 1; this leads to the couplings 



F 2 



(35) 
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All 1PI cactus diagrams contributing to the matrix element in fll8p can be 
generated by the diagrams with full propagators shown in Fig. [3J Each of 
the tadpole loops in Fig. [3] must contribute a factor of N, to be leading in 
N, thus it must be generated by the contraction of a <p T (f) = <j) b (j) b pair, with 
a sum over the flavour index b. All other contractions only give subleading 
powers in N. 

The first term in (|35p then gives only terms with at least one loop integral 
that vanishes since it is odd in momentum. 
We are thus left with the term 



Fd^A* 



yp ~ -Fd^A" 



kbAb 



1 



2F 2 



+ 



(36) 



When the (p b are contracted, they give i 2 times the loop integral A(M 2 hys ) = 

J jT^y I (p 2 — M 2 hys ). Note that the mass in this expression is the physical 
mass and that the propagators in Fig. [3] in the loops are the full propagators. 

We now show that the wave function renormalization vanishes in para- 
metrization 1. The part of the Lagrangian that can produce momentum 
dependence in the full propagator is given by 



kin 



F 2 
2 ' 




(37) 



The corrections to the canonical kinetic term come from the first term, which 
can be rewritten as 

1 (0 a <9^ a )(0 6 <9^ 6 ) 



/'corr 



-"kin 



2F 2 
1 



4> T 4> 

F 2 



(^f)(W)(l + 



+ 



(38) 



In order to have a non-zero loop diagram the derivatives must either both act 
on internal legs or both on external legs. Either way, since the contracted 
legs must have the same flavour, a = b, so there can be no sum over the 
flavour index and thus no leading in N correction. This means that in the 
large iV approximation one has Z = 1 in this parametrization (this will not 
be true in general). 

Putting the terms together, we find the physical decay constant to be 
related to the low energy constants F and M 2 by the simple relation 



phys 



N 
F 2 ' 



l + —A(M 2 hys ) 



(39) 
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To get the LL series one should expand the square root, replace A(M 2 hys ) by 

A = ^ log/i 2 /M 2 hys and express the M p 2 hys in terms of L — j^rps log 

To express this result in terms of the physical -F p h y s and M phys instead, 
we can square f[3T?l) to obtain 

1 + |^(M 2 hys ) = — ^ , (40) 



phys 



2 



which allows us to rewrite ( 13"3l) and (I3"9j) as 



M p 2 hys = M 2 



N — 

\ phys 



phys 



1 - j^A(M 2 hys ) 

phys 1 y 



(41) 



Vacuum Expectation Value 
In a similar fashion one can calculate the leading N vacuum expectation 
value series. Consider the second term in flTj 



F 2 X T $. (42) 
where x T = 2Bq(s° • • • 0). In the first representation this becomes 



2BF^l-^=2BF^l-i^ + ...j. (43) 

In this case there is no need for wave function renormalization. In Fig. H] 
we show how the vacuum condensate is given by the sum of all the tadpole 
diagrams obtained by contracting the C C fields in all possible ways. As 
explained for the decay constant, the leading in N contribution comes from 
the contractions of the same flavour index, i.e. <p c (j) c . Each loop again implies 
the integral A(M 2 hys ). This leads to the following expressions for the vacuum 
expectation value in terms of the low energy constants F, M 2 and in terms 
of the physical -Fphys and M 2 hys 



N — 

V phjs = V Jl + ^A(M 2 hys ), (44) 
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Figure 4: The equation for the vacuum expectation value V^hys- The open thin lined dots 
indicate an insertion of — s . A thick line indicates the full propagator. 



or 

V phys = V ° (45) 

*' ;^(M p 2 hys ) 



phys 



4-3. it re -scattering amplitude 

The meson-meson scattering amplitude denned in ([21]) in the leading iV 
approximation is somewhat more involved. There is ample literature on the 



subject, see for example [15|, |27| and more recently [13|, [l7j, which deal with 
linear or nonlinear massless 0(N + 1)/0(N) sigma models. In the massive 
case there is the additional complication that tadpoles do not vanish. As 
for the previous observables, in order to be leading in N, each momentum 
loop must correspond to a sum over the N flavours. We deal with all the 
generated cactus diagrams in three steps. 

First, we consider all insertions on a meson line that do not carry away 
momentum and/or flavour. They can be dealt with simply by using the full 
propagator obtained earlier. 

Next we deal with effective four-meson couplings, described by the left 
hand side (lhs) of the Eq. in Fig. [51 These can be produced by resumming 
all the generalized tadpoles shown on the rhs in Fig. El just as we did for 
Mp hys , F phys and for Vp hys , see Figs. [21 [3] and H respectively. 

In the massive non-linear sigma model case, the effective four-0 couplings 
come from both the kinetic and the mass part of the Lagrangian. In the first 
representation these may be written as 




£( n>4 )0 = _ v^^/x _ f^+F*M^l - ^-F 2 M 2 + ^0 T 0. 

(46) 
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Figure 5: The effective vertex for meson-meson scattering. The thick lines correspond to 
the full propagator produced by the graphs in Fig. [2j The large dots are vertices from the 
Lagrangian (fT4|) . 

Note that we have removed the kinetic terms that give the lowest order 
propagator. The first term may be expanded into 

4->* = - ± (WW) £ (gf. (47, 

1 F 2 

The loops appearing on the rhs of the equation in Fig [5] may be treated in 
the leading N limit as before. The derivatives cannot both appear in the 
loops at leading order in N for the same reasons valid for wave function 
renormalization. 



The thus leads to an effective vertex 



.(n>4)# 

unun icauo tu an cncuuivc vt 
2F 2 1 + y 



(48) 



where y = -p2A(M 2 hys ). The loop integral A(M 2 hys ) is again produced by the 



b T (b = (b c (b c contractions. 



may be expanded with = E~ =0 ( ' n j (-^) • As 

before, the loops must each come from one (f) T (j) pair, but now we have to 
take into account the number of ways in which the T pairs can be attached 
to the four external legs. For a term (<f) T (f)) n there are n{n — l)/2 ways to 
select four outer fields and to contract the remaining n — 2 pairs. If each 
contraction leads to a factor — A(M 2 hys ), then for each flavour the (0 T 0) n 
term will contribute O T 0) 2 /2 n(n - l)[-A(M 2 hys )/F 2 ] n . The £^| s 4)0 will 

then contribute (0 T 0) 2 /2E^° =O ( ^ 2 J n{n-l)[-N A{M 2 iys )/F 2 ] n , which is 



n 
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Figure 6: The remaining diagrams for A(s,t, u). The vertex is the one obtained from the 
effective vertex (HOj with S ab 6 cd . 



the second derivative of yjl - NA(M 2 hys ) /F 2 with respect to NA(M 2 hys ) / F 2 
The effective four meson vertex coming from the mass term is thus 



M 2 



8F 2 (l + ?/)( 3 / 2 ) 



(49) 



Using fl33l) and fl39|) one may write the total effective vertex which has the 
same form as the lowest order vertex but with physical quantities 



2F phys 



4> a d^4> a 4> h d^ b 



o p2 1 
5 ^phys 



Lb ib 



(50) 



This effective vertex corresponds to the first diagram in Fig. [6] and leads to 
an amplitude (p a +Pb) • (Pc+Pd)—M 2 hys . Each of these vertices is the building 
block for the remaining diagrams in Fig. |6j 

We now concentrate on the A(s, t, u)S ab S cd part of the amplitude defined 
in (1241) . The corrections to (p a +Pb) • (p c + Pd) — ^ph ys are generated by the 
fish-like diagrams given in Fig. [61 This sum is very similar to the sum over 
bubbles used in the NJL model [25( . The structure of the vertex is such that 
for leading iV it does not depend on the loop integral. Each flavour loop in 
Fig. [6] thus adds a factor 



N(s - M 2 hys ] 
IF 2 

phys 



J?(M p 2 hys ,M p 2 hys ,s 



(51) 



with s = (p a + pb) 2 and B the standard two-propagator loop integral 



B(m 2 m 2 ,p 2 ) = - [ 
% J 



d d q 



(2vr) d (q 2 -m 2 )((q-p) 2 



(52) 
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Figure 7: The recursive equation which generates all the meson-meson scattering diagrams. 
The vertex is the effective vertex of ([50]) . The double line is the full results for A(s,t,u). 
The single thick line is the full meson propagator. 

The factor 1/2 is from the symmetry factor in the loop. The sum of diagrams 
forms a geometric series which becomes 

A(s^u) = S —^ — W2 ) -. (53) 

phys 1 - f -r£r B KV' M p 2 hys , s) 

This expression is in agreement^ with both the results in ( 163|) and what was 



found by [13] in the M 2 — > limit. 

An alternative way to resum the diagrams in Fig. [6] is with a recursion 
relation, as depicted in Fig. [71 Let's denote A(s,t,u) by a thick double line 
and a wave function renormalized leg by a thick single line (remember that 
in the large N limit Z = 1). The lhs of the equation is then is result we 
sought after. On the rhs we have the renormalized effective 40 vertex plus 
the A(s,t,u) multiplied the renormalized fish diagram. 

As for the mass case, solving the equation by first writing the lowest order 
expression on the rhs, and reinserting the solution into the equation and so 
on generates the whole set of diagrams in Fig. [6j 

4-4- Vector and Scalar form factors 

The vector and scalar form factors Fy and F$ in the large iV limit are 
calculated in much the same manner. One constructs effective vertices and 
then sums the diagrams. 

The result for the vector form factor is particularly simple. The vector 
form factor couplings come from £ kin , in particular from F 2 (Q T v fl d fJ 'Q — 



5 Compared with [27[ we found an extra factor 1/2 in front of the B(M^ hys , M^ hys , s) 
function coming from the symmetry factor of the loop. Note that they only worked to 
first order in the mass. 
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Figure 8: The fish diagrams for the scalar form factor. 

9 M $ T f M $). The effective vertex for the first parametrization vf [0(9^] is the 
same as the lowest order vertex and because of the antisymmetry in flavour 
indices of v<7 there are also no diagrams similar to those of Fig. [SJ The full 
leading order in N result is thus 

F v (t) = 1 . (54) 

The scalar current comes from £ km , which was discussed earlier for the VEV. 
The sum of all the tadpole diagrams leads to the effective vertex 



Bs*- 



(55) 



In the scalar case however, one must also consider fish-like diagrams, see 
Fig. [HJ The arguments used for resumming these diagrams used for 
scattering still apply. 

The full result for the scalar form factor is thus 



Fs(t) 



F s (0) 



(M p \ ys ,M p \ ys ,tV 



phys 



with 



^(0) 



phys 



F 2 

phys 



(56) 



(57) 



Here we have used the earlier results to write the expression in its simplest 
form. 



5. Leading Logarithmic series for 0(N + 1)/0(N) 

In this section we describe the calculation of the leading logarithms and 
quote results for various physical quantities. For a given observable O p h ys we 
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i 


<2j for N = 3 


a, for general iV 


1 
2 
3 
4 

5 


-1/2 
17/8 
-103/24 
24367/1152 

-8821/144 


1 — 1/2 iV 

7/4 — 7/4 iV + 5/8 N 2 

37/12 - 113/24 N + 15/4 N 2 - N 3 

839/144 - 1601/144 N + 695/48 N 2 - 135/16 N 3 

+231/128 N 4 

33661/2400 - 1151407/43200 N + 197587/4320 N 2 
-12709/300 N 3 + 6271/320 N 4 - 7/2 N 5 



Table 1: The coefficients a, of the leading logarithm L l up to i 
mass [3. 



5 for the physical meson 



can write the leading logarithm expansion in many equivalent ways. The two 
we will use are of the form 



O 



phys 



O {l + a l L + a 2 L 2 + 



(58) 



with L defined in (1291) . We will also expand alternatively in the physical 
quantities 

Ophys = O (l + ClL 



phys 



+ C 2-^phys 



+ ■ 



(59) 



with L phys defined in (!30|) . In both cases we have chosen the mass scale in 
the lo gari thm to be the corresponding mass. 



In 



18] we described how to systematically take into account all the nec- 
essary diagrams for the renormalization of the <ft mass up to five loops. At 
each order new diagrams are necessary. To renormalize the mass at one loop 
for instance, one must consider the C° 4<p vertex and contract two of the legs. 
As the loop order grows, so does the number of outer legs one must consider 
and the number of possible one loop diagrams contributing. We give the 
related discussion here for the decay constant. The actual calculations were 
performed by using FORM 19| extensively. 



18 and 



5.1. Mass 

The coefficients cij of the logs L l for the mass were calculated in 
are here reproduced for completeness in Tab. [TJ The coefficients in terms 
of fully physical quantities can be derived by using the results for the decay 
constant given below. They are given in Tab. |2j The leading logarithm for 
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i 


d for N = 3 


Cj for general iV 


1 
2 
3 
4 

5 


-1/2 
7/8 
211/48 
21547/1152 

179341/2304 


1 - 1/2N 

-1/4 + 3/4 V- 1/8 N 2 
-5/12 + 7/24 iV + 5/8 iV 2 — 1/16 N 3 
347/144 - 587/144 N + 47/24 N 2 + 25 /48 N 3 
-5/128 A^ 4 

-6073/1800 + 32351/2400 N - 59933/4320 N 2 
+224279/43200 N 3 + 761/1920 - 7/256 N 5 



Table 2: The coefficients Cj of the leading logarithm L 
meson mass. 



phys 



up to i = 5 for the physical 



the masses for A" 
two-loop order in 



= 3 at two-loop was first calculated in [3] and later to full 
28l . l29j. Our results agree with those. 



In [18| we noticed that the expansion of M 2 /M 2 hys in -£/M phys converged 
faster than the expansion of M 2 hys /M 2 in terms of L. This was true for both 
the large A" result and the general A" case. From the large A" result in ( 141 p we 
would have naively expected to see a similar improvement in the expansion 
of M 2 hys /M 2 in terms of £phys- Looking at the coefficients of Tab. [2] one can 
see this is not the case. For completeness we also looked at the series of 
M 2 /M 2 hys in terms of L p h ys - The coefficients are of similar size as those in 
Tab. El 

We can now use these results to check the expansions and how fast they 
converge. In 18j the x-axis in Figs. 6(a) and (b) was unfortunately mis- 
labeled. It should have been M [GeV] instead of M 2 [GeV 2 ]. We have 
therefore included a similar figure again. We chose F = 0.090 GeV and 
fi = 0.77 GeV for the plots presented here in Fig. |9J The expansion can 
also be done in the physical quantities and these we show as a function of 
M 2 hys with -Fphys fixed at 0.093 GeV in Fig. [10J Both cases have a similar 
convergence which is fairly slow for masses above about 300 MeV. 

5.2. Decay constant 

The decay constant -Fphys is defined in ( !IH|) . We thus need to evaluate a 
matrix-element with one external axial field and one incoming meson. The di- 
agrams needed for the wave function renormalization were already evaluated 
in the calculation for the mass [18] since we evaluated the inverse propagator 
there. What remains is thus the evaluation of all relevant 1PI diagrams with 
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(a) (b) 

Figure 9: The expansions of the leading logarithms order by order for F = 0.090 GeV, 
H = 0.77 GeV and N = 3. (a) M 2 hy JM 2 , expansion in L. (b) M 2 /M 2 hys , expansion in 

L M phyB ■ 
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(a) (b) 

Figure 10: The expansions of the leading logarithms order by order for -Fphys = 0.093 GeV, 
H = 0.77 GeV and N = 3. (a) M 2 hys /M 2 , expansion in L phys . (b) M 2 /M 2 hys , expansion 

in -C'phys- 
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(a) (b) (c) 

Figure 11: The diagrams needed up to order 3 for the one particle irreducible diagrams 
with an external meson and an axial current. The wiggly line indicates the axial vector. 
Vertices of order ft 1 are indicated with Q]. (a) The diagram needed at order ft. (b) The 2 
diagrams needed at order ft 2 , (c) The 4 diagrams needed at order ft 3 . 




Figure 12: The diagrams needed for the divergence of the 3- meson 1-axial- vector vertex, 
(a) The 2 diagrams to order ft. (b) the 5 diagrams to order ft 2 

an external a". At order % there is only one diagram, at U 2 there are 2 and at 
order ft 3 there are 4. These are shown in Fig. [TTJ We have not shown them 
but at order % A there are 7 and at h 5 there are 13 diagrams to be calculated. 

To order % it is sufficient to know the lowest-order Lagrangian, but at 
order Ti 2 we need to know the (divergent part of the) vertices coming from 
the Lagrangian of order K with one and three external meson legs and one 
a a . The diagram of Fig. UTT a) gives the divergence of the vertex with one 
meson and one axial vector leg but we also need to calculate the divergence 
of the vertex with three meson legs and one axial vector. This requires the 
diagrams shown in Fig. [121(a). 

To order h 3 , we need still more vertices, we need the divergence of the 
one meson one axial vector leg vertex to order % . These diagrams we have 
already calculated, but we also need the four-leg vertex to order h 2 which 
can be calculated from the diagrams in Fig. [12T b). Inspection of the vertices 
there shows we already have all we need but for the five meson one axial 
vector leg vertex at order %. To obtain that we also need to evaluate all 
diagrams shown in Fi g. 1131 By now, the pattern should be clear, to get the 
mass at order h n in Il8fl we needed the 2 and four-meson vertex at order 
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Figure 13: The 3 diagrams needed for the divergence of the 6-meson vector to order ft. 



h n -\ the 2, 4 and 6 -meson vertex at order Ti n 2 and so on. Here we need for 
the decay constant to order h n in addition the one meson one axial vector 
and three-meson-one-axial- vector vertex at order h n ~ l , the 1, 3 and 5-meson 
plus one axial vector vertex at order h n ~ and so on. Similarly one can see 
that to get the decay constant at order h n , we need to calculate one-loop 
diagrams with up to n vertices. The extension to order h 5 shows that we 
we need to calculate the 18, 26, 33, 26 and 13 diagrams at orders h , . . . , h 5 
respectively for the mass/wave function renormalization and an additional 
27, 45, 51, 33 and 13 diagrams at orders h , . . . , h 5 for the 1PI diagrams with 
an axial-vector. 

We have organized this calculation by first expanding the lowest-order La- 
grangian to the order needed, up to vertices with 12 mesons or 11 mesons and 
one axial vertex. With these vertices we then calculate all 1PI diagrams with 
up to 10 external legs. The divergent part of all needed integrals can be cal- 



culated relatively easily using the technique described in App. Appendix A 



of [18l | . At this stage, the dependence on external momenta is also put back 
as derivatives on the external legs and everything assembled to give the di- 
vergent part at order h for all the vertices with up to 10 meson legs or nine 
mesons plus one axial- vector using (fl2"l) . So we have assembled everything 
we need to calculate the one- loop divergences to order h 2 . The 26+45 dia- 
grams are evaluated and we obtain the divergences at order h 2 using (fT2"]) . 
The process is then repeated up to order h 5 . All of the above steps have 
been programmed in FORM. The CPU time needed increases rapidly with 
the order n one wishes to reach. The Lagrangians at higher orders tend to 
contain very many terms and constructing the diagrams with many external 
legs at higher orders is also extremely time consuming. The CPU time used 
on a typical PC for the mass-divergence to order h n was approximately 0.1 
seconds for h, 0.3 seconds for h 2 , 11 seconds for h 3 , 700 seconds for h 4 and 
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i 


<2j tor iv = o 


Oj for general A~ 


1 
2 
3 
4 

5 


1 

-5/4 
83/24 
-3013/288 

2060147/51840 


-1/2+ 1/2 JV 

-1/2 + 7/8 JV- 3/8 iV 2 

-7/24 + 21/16 iV- 73/48 N 2 + 1/2 N 3 

47/576 + 1345/864 N - 14077/3456 N 2 

+625/192 N 3 - 105/128 N A 

-23087/64800 + 459413/172800 N 

-189875/20736 N 2 + 546941/43200 N 3 

-1169/160 N 4 + 3/2 N 5 



Table 3: The coefficients ai of the leading logarithm L l for the decay constant F p h ys in 
the case N = 3 and in the generic N case. 

30000 seconds for h 5 plus a similar amount for the extra diagrams needed for 
the decay constant. 

We now give the leading logarithms for the decay constant as a function 
of F and M 2 and of the physical -F p hys 

and M 2 hys 
F phys = F (l + axL + a 2 L 2 + ■ ■ ■) , 

F phys = F(l + Cl L phys + c 2 L 2 phys + ■■■) . (60) 

The first five a, coefficients are listed in table [3] for the generic A^ and for 
the interesting case A^ = 3. The equivalent results for the first five q are in 
table HI Note that once the expression of -Fphys as a function of F is known 
one may express the remaining observables as a function of the physical M 2 hys 
and Fphys- This has already been used to calculate the q coefficients in tables 
|5] and H] from the corresponding a*. 

We have plotted in Fig. Q3] the expansion in terms of the unrenormalized 
quantities and in terms of the physical quantities. In both cases we get a 
good convergence but it is excellent for the expansion in physical quantities. 

5.3. Vacuum expectation value 

The expression for the leading logarithms of the vacuum expectation value 
Vphys follows from the definition (fT9~j) . The diagrams needed and the principle 
of the calculation can be derived in the same way as we did for the decay 
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i 


d for N = 3 


Cj for general iV 


1 
2 
3 
4 

5 


1 

5/4 
13/12 
-577/288 

-14137/810 


-1/2+ 1/2 iV 

1/2 - 7/8 N + 3/8 N 2 

-1/24 + 13/16 N - 13/12 iV 2 + 5/16 N 3 

-913/576 + 2155/864 iV - 361/3456 N 2 - 69/64 iV 3 

+35/128 N 4 

535901 /129600 - 2279287/172800 N + 273721 /20736 N 2 
-11559/3200 N 3 - 997/1280 N A + 63/256 N 5 



Table 4: The coefficients c, of the leading logarithm L l phys for the decay constant -Fphys in 
the case N = 3 and in the generic N case. 
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n 


a n for N = 3 


a n for general iV 


1 
2 
3 
4 
5 


3/2 
-9/8 
9/2 
-1285/128 
46 


+1/2N 

+3/4 3/8 N 2 
+3/2iV-3/2A^ 2 + l/2 N 3 

+145/48 N — 55/12 iV 2 + 105/32 iV 3 - 105/128 N 4 
+3007/480 N - 1471 /120 N 2 + 557/40 N 3 
-1191/160A^ 4 + 3/2 N 5 



Table 5: The coefficients at of the leading logarithm L % for the VEV V p hy S in the case 
N = 3 and in the generic N case. Note that the coefficients in front of the first subleading 
N power are often large. 



n 


c n for A" = 3 


c n for general A" 


1 

2 
3 
4 
5 


3/2 
21/8 
75/16 
1023/128 
2669/256 


1/2 N 

-1/4 N + 3/8 N 2 

1/4JV-1/2 N 2 + 5/lQN 3 

3/16 N + 5/24 N 2 - 59/96 N 3 + 35/128 N 4 

-4153/2880 N + 12299/4320 N 2 - 142/135 N 3 

-167/320 N 4 + 63/256 A" 5 



Table 6: The coefficients a of the leading logarithm L l phys for the VEV V p h ys in the case 
V = 3 and in the generic V case. 

constant in Sect. 15.21 The first five coefficients aj defined by 

V phys = -2BF 2 (l + a x L + a 2 L 2 + ••■),, (61) 

are given in table [5] for generic A" and for the interesting case N = 3. The Cj 
for the leading logarithms in terms of physical quantities are given in table [6j 

We have plotted in Fig. [TS] the expansion in terms of the unrenormalized 
quantities and in terms of the physical quantities. In both cases we get a 
good convergence but it is excellent for the expansion in physical quantities. 

5.4- 00- scattering amplitude 

Elastic 7T7r-scattering is a the ideal interaction to test the convergence of 
the ChPT expansion since it only involves the SU(2) L x SU(2) R /SU(2) V 
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(a) (b) 

Figure 15: The expansions of the leading logarithms order by order for fi — 0.77 GeV and 
N = 3. (a) Vphy S /Vb in terms of M 2 , expansion in L with F = 0.090 GeV. (b) V p hy S /Vo 
in terms of M^ h expansion in L p h ys with F p h ys — 0.093 MeV fixed. 



Goldstone bosons. It is also the simplest purely strong interaction process. 
This interaction has indeed been studied precisely for these reasons first by 
24 1 at tree level and then by 0, [si at one- loop level. One would expect 



the series to converge rather rapidly, instead the authors of |3[ found that 
some one-loop corrections were rather large. Specifically, as mentioned in 
the introduction, for the a® scattering length. On the other hand, the one- 
loop LL correction to a\ arises only through F n renormalization [7|, which 
means that the chiral logs are not the main part of the one-loop correction. 
The issue of the convergence of the perturbative ChPT expansion for 00 
scattering is delicate. 

We have above obtained the all-order leading in large N result. From 
the calculations for the mass we have already obtained all the needed diver- 
gences to get the LL to four-loop order for meson scattering. The result can 
be written as expected in the form of ff24|) . This is the first result where 
the choice of scale in the logarithm is not unique. We only quote here the 
expansion in terms of physical quantities and define 



M phys 



A 



u 



-^phys 
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Lm = ia§f; log ^ (62) 

with a generic scale M. For scattering lengths an obvious choice is M = 
Mphy S but in the massless case the choice is M 2 = s. Our result for general 
N is 



F 2 

phys 

Mi 



2 



A(s, t, u) = 5 - 1 + L M [(1/6 A 2 - s 2 + 1/2 N § 2 ) + (11/35 — N s) 



phys 

+ (-8/3 + 1/2 iV)j + L 2 ^ [(5/96 5 A 2 + 181/288 S 3 + 5/96 N 5 A 2 
- 163 /288 N 5 3 + 1/4 iV 2 § 3 ) + (-5/12 A 2 - 91/36 S 2 + 5/12 N A 2 
+29/18 iV s 2 - 3/4 N 2 5 2 ) + (25/18 s + 23/18 N s + 3/4 iV 2 5) + (4/3 
-29/12 iV- 1/4 N 2 )] +L 3 M [(361/17280 A 4 - 317/12960 S 2 A 2 

-21319/51840 s 4 - 203/17280 iV A 4 + 229/6480 N § 2 A 2 
+28081 /51840 A^s 4 + 1 /160 N 2 A 4 + 11/1440 N 2 5 2 A 2 - 33/80 N 2 § 4 
+ 1/8 N 3 S 4 ) + (-1901/25920 5 A 2 + 51869/25920 S 3 + 3073/25920 A^s A 2 
-49573/25920 A^s 3 + 41 /288 N 2 s A 2 + 8467/4320 A^ 2 s 3 — 1 /2 N 3 s 3 ) 
+ (1283/6480 A 2 - 907/720 S 2 - 2503/2160 N A 2 - 7193/6480 A r s 2 
+43/60 N 2 A 2 - 3257/1080 N 2 s 2 + 3/4 N 3 s 2 ) + (-1189/1620 5 
+2111/810 AT s + 607/108 A^ 2 s — 1/2 A^ 3 s) + (17/810 + 457/180 N 
-22/5 N 2 + 1/8 N 3 )] + L A M [(1451/1244160 5 A 4 + 64 57/103 6 8 5 3 A 2 

+61781/248832 s 5 + 143893/12441600 A^s A 4 - 77957/2073600 A" s 3 A 2 
-5387831/12441600 A^s 5 — 9089/1382400 N 2 s A A + 5531 /230400 N 2 5 3 A 2 
+5592583/12441600 N 2 § 5 + 1/256 N 3 § A 4 - 1/3840 N 3 § 3 A 2 
-5267/21600 N 3 § 5 + 1/16 N 4 § 5 ) + (-6493/77760 A 4 
+9023/103680 S 2 A 2 - 684899/518400 S 4 + 43523/345600 N A 4 
-203777/1036800 N § 2 A 2 + 20749/12150 A^s 4 - 19091/345600 N 2 A 4 
+146869/1036800 N 2 s 2 A 2 - 1840297/777600 N 2 § 4 + 7/320 N 3 A 4 
+143/5760 N 3 s 2 A 2 + 110897/86400 N 3 § 4 - 5/16 N 4 5 4 ) 
+ (-680609/1555200 5 A 2 + 23719/103680 S 3 - 331117/1555200 A^5 A 2 
+2894959/1555200 A^s 3 + 16621 /86400 A^ 2 s A 2 + 2812631/777600 N 2 s 3 
+77/288 A^ 3 5 A 2 - 153377/86400 § 3 + 5/8 N 4 5 3 ) + (39629/15552 A 2 
+88013/129600 S 2 - 186451/129600 N A 2 - 272671/77760 A^5 2 
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Uq 


a 2 




b 2 


d 1 

Cl 

c 2 
c 3 

C 4 




-2M\ 

pnys 


8M\ 

pnys 


-AM 2 , 

pnys 


32ttF' 2 . 

phys 

9/2 
857/42 
153211/1512 
41581/84 


32vrF 2 , 

phys 

-3/2 
-31/6 
-7103/216 
-7802/45 


32wF 2 , 

pnys 

26/3 
1871/36 
2822/9 
744463/43 


32wF' 2 , 

phys 

-10/3 
-169/36 
-352/9 
-1309703/6480 



Table 7: The coefficients Cj of the leading logarithm series dj phys = d\ tree {\ 



4 



7TTT — >7TTT 



6q, b\ in the case N 



C l^phys + 

for the scattering lengths, d\ = Oq, dp, and for the slopes 
3. All in units of M p h ys - 



-9227/5400 A^ 2 A 2 - 48067/6075 A^ 2 s 2 + 131/120 N 3 A 2 
+2017/3600 A^ 3 § 2 - 5/8 A^ 4 5 2 ) + (667007/48600 5 
-1109347/129600 Ns + 369719/43200 N 2 s + 2467/432 A^ 3 S 
+5/16 iV 4 5) + (-12349/864 + 102659/10800 + 36097/10800 N 2 
-2887/480 N 3 - 1/16 iV 4 )" . 



(63) 



As for all other quantities we see large subleading in N corrections. 

From the result (16"3"|) we can obtain the different amplitudes T 1 defined in 
(I2"6"j) and project on the partial waves using (12"7|) . The scattering lengths and 
slopes as defined in fl28l) can then obtained as well and we get the LL L l phys 
in terms of the physical M 2 hys and F p h ys . We have LL contributions for all 
a\ up to I = 5 and to the slopes up to £ = 4. These we have all calculated 
for general N. We give the expansion in 



do „(,„,„ = d t tree (1 + ci££ 



L i,phys 



phys 



+ C 2-^phys + 



(64) 



for the S- wave scattering lengths and slopes and only quote the phenomeno- 
logically relevant case of N = 3. As mentioned above a clear choice for the 
arbitrary scale in the logarithm is here the physical mass. The lowest or- 
der result and the expansion coefficients are given in table [71 The one- and 
two-loop results agree with the earlier published ones @, 0,130 



6 The two-loop coefficients agree with those of Q except for b^. We have checked that 
using the full result from [29[ and (3.13) in reproduces our result. 



31 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 
Mp h>s [GeV z ] 



(a) 
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Figure 16: The expansions of the leading logarithms order by order for // = 0.77 GeV and 
N = 3. (a) a$/a$ troc (b) a§/a§ troe in terms of M% hys , expansion in L p h ys with F ph ys = 
0.093 MeV fixed. 

We have plotted in Fig. [16] the expansion in terms of the physical quanti- 
ties of <2q and Oq. There is an excellent convergence for mass up to 0.2 GeV 
but above 0.3 GeV it becomes rather slow for a[j. For a 2 , it is somewhat 
better but also rather slow at the higher masses. 

In Fig. [17] we plotted the expansion in terms of the physical quantities 
of the slopes 6° an d b$. There is an excellent convergence for mass up to 
0.25 GeV but above the convergence slower for b$. b^ converges better but 
also rather slow at the high mass end. 

In the massless case we can obtain the coefficients at higher orders also 
with different methods [HI, and our result agrees with those. We can also 
use our result in this limit to test the often used elastic unitarity arguments. 
The partial waves T/ satisfy under the assumption of elastic unitarity 



In the massless case the scale of the logarithm should be related to s and we 
know that this should come in the combination /(/i 2 /s) = In |/i 2 /s| + iir8(s). 
The leading logarithm part can be written as 




rpl 2 



(65) 



(66) 
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(a) (b) 

Figure 17: The expansions of the leading logarithms order by order for \i = 0.77 GeV 
and N — 3. (a) &o/&o trco ^) ^o/^otrcc m terms of M£ hys , expansion in L p h ys with F p h ys = 
0.093 MeV fixed. 



In the chiral limit the interaction must vanish at s = 0. Inserting f l66l) in 
(|65|) and keeping only the leading logarithm part allows to determine all the 
higher coefficients in terms of the lowest order one. 




< = ■ (67) 



Note that in the massless case there is really no regime where elastic unitarity 
( )65|) is valid. We are simply testing here how much of the leading logarithm 
in this case follows from the so-called righthand two-body cut. The result 
( 16"T|) can be written explicitly 

T/ = & . (68) 

1 - ^/s) 

In table [8] we compare for N = 3 the exact coefficients derived from (1631) with 
those from ( |68l) for the expansion 



' (69) 

It is clear from the table that the assumption of elastic unitarity does not 
give a good approximation to the LL in the chiral limit. 
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n 


T 2 


el 


i 


i 


-l 




WttF 2 , 

phys 


32wF 2 , 

phys 


32vrF 2 , 

phys 




exact elastic 


exact elastic 


exact elastic 


fi 

h 
h 
U 


25/18 1 
25/18 1 
18637/9720 1 
540707/291600 1 


1/2 
5/4 1/4 
-901/3240 1/8 
207871/170100 1/16 


-10/9 -1/2 
35/36 1/4 
-50077/48600 -1/8 
134077/145800 1/16 



Table 8: The coefficients e\ and /, defined in for the exact LL in the chiral limit and 
those derived using the assumption of elastic unitarity for the case N = 3 and T§, T-j 1 and 



T 2 
1 o ■ 



5. 5. Form factors F$ and Fy 

The vector and scalar form factors were defined in [20] and [22j The pro- 
cedure to find the LL for this observable follow the same lines of the one for 
the decay constant, with the difference that in representation 1 the vertex 
between the meson fields and the vector current is simply given by 

£ int = (1/2) vf [Vfafa - fadvfo] . (70) 

For the wave function renormalization one may again use the results obtained 
for the mass calculation. As in the previous subsection we express here the 
results in terms of t — t/M 2 hys and Lm with a scale Ai 2 some combination 
of t and Mp hys . The result for Fy to four- loop-order for the LL is: 

F v (t) = l + L M \l/Qt\ + L 2 M \t (-11/12 + 5/12 N) + P (5/36 - 1/24 N) 



M 



t (+1387/648 - 845/324 N + 7/9 N ) + t 2 (-4007/6480 



+3521/6480 N - 29/180 N 2 ) + P (+721/12960 - 47/1440 + 1/80 N 2 



+L 



M 



t (—44249/15552 + 222085/31104 — 55063/10368 iV 2 

+127/96 A^ 3 ) + t 2 (+349403/155520 - 15139/4860 N + 86719/51840 A^ 2 
-199/480 A^ 3 ) + t 3 (-85141/155520 + 885319/1555200 N 
-5303/19200 A^ 2 + 21/320 A^ 3 ) + t 4 (+4429/103680 - 57451/1555200 N 
+289/14400 A^ 2 - 1/240 N 3 )] . (71) 

Note that Fy(0) = 1 as it should be. We can extract from this the expansion 
for the radius and curvature defined in ( |23i) . These are given in table [9] in 
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n 


(r 2 )v 


cv 


(r 2 >s 


cs 


Cl 


1 




6 




c 2 


2 


1/72 


-29/3 


43/36 


c 3 


853/108 


-71/162 


-581/54 


-727/324 


c 4 


50513/1296 


-25169/7776 


-75301/648 


4369/810 



Table 9: The coefficients of the expansion in L p h ys in the expansion of the radii (r 2 )v,s 
and the curvature cy,s in the iV = 3 case. 



terms of the expansion in L p h ys for the physical case N = 3. The general 
coefficients can be easily derived from (J7TJ). The dash indicates that this 
cannot appear to a given order for the LL. The result up to two-loop order 



30 . We 



agrees with the LL extracted from the full two-loop calculation 
do not present numerical results for the vector form factor since these are 
dominated in the physical case N = 3 by the large higher order coefficient 
contributions, see e.g. [2I. [3o(] . 

For the scalar form factor Fg(t) defined in ( 1201) we have already done the 
calculations we need to four-loop-order during the calculation for the VEV 
Vphys to five-loop-order. We write the result in the form 



F s (t) = F s (0)f s (t) 



(72) 



The value for Fs(0) can be obtained from the calculation or via the Feynman- 
Hellman theorem in (|21|) . The latter allows for the <Zj coefficients for Fg(0) 
expanded in terms of L to be derived easily from table [TJ We have checked 
that both methods agree. In table [10] we quote the coefficients Cj of 



F s (0) = 2B[1 + c x L + c 2 L 2 + 



(73) 



for the case N = 3 and general N. 

We can now express fs{t) using the same notation as for Fy(t). 

f s (t) = 1 + L M \i (-1/2 + 1/2 N)] + L 2 M \t (1/18 + 7/36 N- 1/4 iV : 



+t 2 (1 1/72 - 29/72 N + 1/4 N 2 



T 3 



t (599/648 - 181/324 N 



-53/108 N 2 + 1/8 iV 3 ) + t 2 (275/1296 - 427/648 N + 301 /432 A^ 2 
— 1/4 A^ 3 ) + t 3 (-89/864 + 7/24 N - 271 /864 A^ 2 + 1/8 N 3 ) 



+L 



M 



t (-10981/3888 + 37373/7776 N - 3733/2592 A^ 2 - 23/48 A^ 3 
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Cj for N = 3 


Cj for general iV 


Cl 

c 2 
c 3 
c 4 


-1 
31/8 
65/6 
76307/1152 


2-N 

5/A-l/AN + 3/8 iV 2 

-5/3 + 41/12 iV + A^ 2 - 1/4N 3 

655/144 - 901 /144 N + 341 /48 iV 2 + 17/48 N 3 

+11/128 A^ 4 



Table 10: The coefficients Cj of the leading logarithm L l phye for the scalar form factor at 
zero momentum transfer F$(0) in the case N = 3 and in the generic N case. 



-1/16N 4 



-22859/28800 + 89951/48600 N - 823067/777600 N 2 



-4021 /21600 N 3 + 3/16 N 4 ) + P (-959/32400 + 82529/259200 N 
-1421 /2025 N 2 + 51877/86400 N 3 - 3/16 N 4 ) + £ 4 (76459/1555200 
-70997/388800 N + 423961/1555200 N 2 - 121 /600 N 3 + 1 /16 N 4 ) . 

(74) 

From (lT4"j) we can derive the expression for the scalar radius and curvature 
defined in ff23l) . The expansion coefficients Cj in terms of the physical loga- 
rithm Lphy S are given in table [9] for the physical case N = 3. The general 
case can be easily obtained from (1741) . The coefficients of the LL extracted 
from the full two-loop calculation of 30[] agree. 

In Fig. [18] we plotted the expansion in terms of the physical quantities of 
the radius and curvature. There is an excellent convergence for masses up to 
0.3 GeV but above the convergence is slower in both cases. 



6. Conclusions 

In this work we extended our previous work on the massive nonlinear 
0(N + 1)/0(N) sigma model to many more observables. We calculated the 
leading logarithmic corrections to the decay constant and the vacuum ex- 
pectation value to five-loop-order and to the scalar and vector form factors 
and meson-meson scattering to four-loop order for generic N. We used these 
results to extract scattering lengths and slopes and have presented the phys- 
ically most relevant cases for A^ = 3 of these. Results for all other cases have 
been obtained but would have added significantly to the length of the paper. 

Our original hope had been to find a pattern behind the coefficients of 
the LL and to make an all order conjecture for the leading LL. We succeeded 
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Figure 18: The expansions of the leading logarithms order by order for \i = 0.77 GeV 
and N — 3. (a) M^ hys (r 2 )s (b) Afp hys cs in terms of M^ hys , expansion in L p h ys with 
F phys = 0.093 MeV fixed. 



in deriving such an expression for the leading term in the large N limit but 
we found no general expression. 



The large N approximation, as we already noted in [18| , is a surprisingly 



poor approximation of the LL series for all of the observables we considered. 
For example, looking at the five-loop results, the first neglected term, the 
iV 4 term, often has a large coefficient compared with the N 5 term. For this 
term to be negligible, i.e. a 10% correction of the leading term, N must be 
large, in some cases N > 20. This is understandable if one considers that 
the subleading N l<n terms in the coefficients come from non-cactus diagrams 
and different flavour routings of the cactus diagrams. Though each of these 
diagrams is suppressed by 1/N in the large N limit, the number of diagrams 
and the number of ways to route the flavour structure seem to produce large 
coefficients for the subleading in N terms. 

We have also performed some numerical test of the convergence with pa- 
rameter values of the range needed for two-flavour ChPT. For masses around 
0.1 GeV the convergence for all quantities studied is excellent. It is reasonable 
for most quantities up to about 0.3 GeV as can be seen on the various plots 
but becomes unstable around 0.4 GeV for some of the quantities studied. 
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